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By  an  extension  of  ordinary  geometrical  optics 
(or  acoustics)  the  intensity  of  the  reflected  and  trans- 
mitted fields  due  to  a  point  source  in  the  presence  of 
an  arbitrary  interface  between  two  media  is  found.  A 
particular  consequence  of  the  solution  is  the  general 
lens  and  mirror  law  and  the  equations  for  the  caustic 
surfaces. 


1.  Introduction 

The  calculation  of  the  electromagnetic  or  acoustic  field  at  any 
point  in  space  when  two  different  media  are  present  necessitates  a  solu- 
tion of  Maxwell's  equations  or  the  wave  equation  of  acoustics  with  appro- 
priate "boundary  conditions  in  each  case.  Exact  solutions  have  teen  ob- 
tained  only  for  the  simplest  configurations  '  and  approximate  solutions 
have  "been  found  by  special  means  in  some  other  cases  (e.g.  when  one  of  the 
media  is  an  arbitrarily  curved  thin  shell^'  ;. 

A  general  procedure  for  obtaining  an  approximate  solution  to  an 
electromagnetic  or  acoustic  problem  is  the  method  of  geometrical  ox>tlcs 
or  acoustics.  By  this  it  is  meant  that  the  field  quantities  propagate 
along  rays  which  are  determined  by  the  Fermat  principle  of  least  time, 
and  that  these  rays  obey  the  laws  of  reflection  and  refraction  at  the 
interface  between  two  different  media.  Uurthermore,  in  this  method  it 
is  also  possible  to  determine  the  field  components  themselves"5*  -  some- 
thing which  has  not  ordinarily  been  done  in  geometrical  optics  or  acoustics. 
Thus,  for  example,  the  reflected  and  transmitted  field  components  at  an 
interface  are  related  to  the  incident  field  components  by  the  well  known 
Fresnel  formulae  or  the  correspond!  r.g  formulae  in  acoustics.  Also  the 
field  components  vary  along  a  ray  inversely  as  the  square  root  of  the 
velocity  and  of  the  area  of  the  normal  cross  section  of  an  infinitesimal 
tube  of  rays  containing  the  ray  in  question.  These  results  enable  one  to 
determine  the  geometrical  optics  or  acoustics  field  completely.  It  is  to 
be  expected  that  this  solution  will  be  an  adequate  approximation  to  the 
full  solution  only  at  very  high  frequencies  (i.e.  at  wavelengths  small 
compared  to  the  dimensions  of  the  problem).* 

*  As  a  matter  of  fact  R.  K.  Luneberg  has  shown  quite  generally  that  the 
first  terms  in  asymptotic  expansions  of  the  space  dependence  of  the 
electric  and  magnetic  field  vectors  is  that  given  by  geometrical  optics. 
This  material  will  appear  in  a  forthcoming  report. 

""Phase  is  determined  by  the  path  length  along  a  ray. 


In  the  present  Investigation  an  incident  field  due  to  a  ooint 
source  is  assumed  to  impinge  upon  an  arbitrarily  curved  or  flat  interface 
which  separates  two  homogeneous  and  isotropic  media.  The  magnitudes  of 
the  reflected  and  transmitted  field  components  at  any  point  are  deter- 
mined on  the  "basis  of  geometrical  optics,  i.e.  "by  applying  the  theory 
outline  above.  For  the  case  of  reflection  the  problem  has  already  been 
solved-'9',  but  it  is  included  here  because  this  can  be  done  with  no  addi- 
tional difficulty.  The  transmission  problem  has  been  treated  by  the 
Kirchhoff  method^*  and  tho  present  calculation  constitutes  a  check  on 
that  solution.  The  check  for  the  Kirchhoff  method  applied  to  reflection 
is  contained  here  as  well  as  in  reference  3* 

2.  Formulation  of  the  Proble: 

A  point  source  is  assumed  to  be  located  at  a  point  (xi,yi,zi) 
in  a  homogeneous  isotropic  medium  with  propagation  speed  VT .  The  surface 
z  =  z(x,y)  separates  this  medium  from  a  different  homogeneous  isotropic 
medium  with  propagation  spaed  Vm« 

The  magnitudes  of  the  reflected  and  transmitted  field  components 
at  the  surface  can  be  found  in  terms  of  the  incident  field  components  and 
the  Fresnel  formulae  (only  for  those  surface  points  which  can  be  connected 
to  the  source  by  a  straight  line  segment  lying  in  the  first  medium).  Let 
E(x,y,z)  represent  the  amplitude  of  any  field  component  at  the  point 

(x,y,z)  due  to  reflection  or  transmission  from  the  surface.  Then  from 
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geometrical  optics  we  have  the  relation 


EU.y.z)  =  SjU'.y'.a').^ 


do- 


where  (x',y',s')  is  the  point  in  which  a  reflected  or  transmitted  ray 
through  (x,y,z)  intersects  the  surface  z»  =  zCx'.y').  S.Cx'.y'.z')  is 
the  reflected  or  transmitted  field  component  at  the  surface,  where  i  =  1 
for  reflection  and  i  =  2  for  transmission.  These  components  can  be  cal- 
culated, as  previously  stated,  from  the  incident  field  component  at  the 


surface.   The  quantity  dcr  is  the  area  in  which  an  infinitesimal  tube  of 
reflected  or  transmitted  rays  containing  the  ray  through  (x,y,z)  and 
(x'.y'.z1)  cuts  a  plane  p  normal  to  this  ray  at  the  point  (x,y,z),  (see 
Pig.  2).  d  cris  the  area  enclosed  by  this  same  tube  of  rays  on  a  plane 
p!  normal  to  the  ray  in  question  at  the  point  (x'.y'.z1).   In  the  limit 

of  infinitesimal  areas  the  ratio  dC/d<S*i3  just  the  Jacoblan  J(2— )  0f 

P  * 
the  transformation  established  by  reflected  or  transmitted  rays,  which 

maps  the  plane  p1  into  the  plane  p.  Hence  we  may  write 


(1) 


E(x,y,z)  =  W2  (21)  E«(x',y',: 


3.  Calculation  of  the  Jacobia] 

Let  E(x',y{,z')  be  a  unit  vector  normal  to  z'  =  z(x',y»)  and 
pointing  into  the  first  medium;   l(x,,yl ,z')  is  the  unit  vector  pointing 
from  (x'.y'.z1)  to  the  source  at  (x^y-j^z-^).  The  direction  of  the  re- 
flected ray  at  (x'.y'.z1)  is  given  by  the  unit  vector  'iCx'.y'.z')  and 
that  of  the  transmitted  ray  by  T(x8,y!,z')  where,  from  the  laws  of  re- 
flection and  refraction  [see  Appendix J  : 

(2)  ^U'*'..')  =  -  n1-1  I  -  p-1^)  -(-l)1{l^1-1>[^(!*^/2]t 

(i  =  1.2) 

sin#  i 
Here  n  -     sin^         and  i  =  1  for  reflection,    i  =  2  for  refraction  or  trans- 
mission.    The  angles  and  unit  vectors   are  shown  in  Fig.   1. 


(*.,»■/.)• 


Pig.   1 


(3) 


— > 
The  unit  vectors  I  and  U  have  the  components: 

*i  -  r-s       y,  -  y'  z.   -  z» 

l(x«  yi  zO  =  / ,-= ,—        1 1 

'.ACx'.y'.a')   '  ^U>,y',z')  '  ^x'.y'.z'), 

-  5z'       -  gz' 

"f(x!  y»  -»)  =  (      5x8  gy' 

VJKxi.yt..!)  '^(x.,ytz.)  '  ^(x.,y.tZ,)  I* 

where: 

J2(x',y',z«)  =  (sfx«)2  +  (y^-y')  ♦  (^  -  z«)2 
and 

The  coordinate  system  must  be  so  chosen  that  N  points  into  the  first 
medium. 

For  an  arbitrary  point  (x'.y'jZ*)  on  the  surface  the  equation  of 
the  reflected  or  transmitted  ray  through  this  point  is 

x  -  x»  „  y  -  y» 


*?  •_  r? 


i  i 

x  y         z 

and  thus: 

**  (x».,y»,«») 


x  =  x'  +  (a  -  a*) 


*T  (x',y\s«) 

(U) 

V_  (x«,y«,s«) 

y  =  y'  f  (z  -  z>)  -5-i 

1Tz(x«,y«,z') 

where  the  subscripts  indicate  components  of  the  vector  T  . 

For  a  fixed  value  of  z  equations  (U)  are  the  mapping  of  the  sur- 
face S  given  by  z »  =  z(x',y')  onto  the  plane  TT,  given  by  z  -  const  (see 
Fig.  2).  From  these  equations  we  can  compute  the  Jacobian  of  this  trans- 
formation, namely, 

(5)      -v    J    (?     I973T       9F  a^-J 


^tf.,2.) 


Fig.  2  (Transformations  for  Transmission) 


The  Jacohian  which  occurs  in  equation  (l)  can  now  he  calculated  since 


'<£■>  ■  «P  ^  J(f  >  J'-P 


p     "  'Tf'  "         -S     •         TT         -p- 

where  TT1    is  a  plane  through  (x'.y'.z')  and  parallel  to  tt,  p  and  p'  are 

the  normal  planes  previously  defined  and  the  notation  represents   the 
Jacohians   of  the  transformations  between  indicated  surfaces   determined 
"by  means  of  the  reflected  or  transmitted  rays.     However,J(i-y)=  J     (— ) 

since  the  planes  are  parallel  "by  pairs  and  thus  corresponding  angles  used 
in  defining  the  transformations  are   equal.   Thus  we  have: 

j(El)  =  j(Ili)  j(l) 

p  S  TT 

How  place  the  origin  of  coordinates  at   the  point    (x^y'.z1)  with 

the   z-axis  normal  to  S  and  positive   into  the  first     medium.     Further  let 
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the  x-  and  y-axes  he  parallel  to  the   directions  of  principal  curvature 

at  this  point.     The  plane  Tt'   now  "becomes  the  xy  plane  (z=C)   and  is  the 
tangent  plane  to  S  at   the  new  origin.     The  equation  of  the  surface  S  ex- 
panded into  a  Taylor  series  around  the  origin  "becomes^ 


(6) 


:«  =  a(x»)     +  h(y')       + 


where  2a  and  2b  are  the  principal  curvatures  of  the  surface  at  the  origii 
and  the  primes  denote  the  coordinates  of  a  surface  point  relative  to  the 
new  coordinate  system. 

TTS  \ 

In  these  coordinates  J(— )  is  the  Jacobian  of  the  mapping,  by 
reflected  or  transmitted  rays,   of  the  surface  on  the  xy  plane.      This 
mapping  is  given  by 


(7) 


x  ■  x»  -  «'  tan    6(x«,y') 
y  =  y*   -  z»   tan  (b  (x',y' ) 


where     ©(x',y')  and  <£  (x'.y1)  are  the  angles  between  a  ray  reflected  (or 
transmitted)  at   (x'.y'.z1)  and  the  x-  and  y-axes,    respectively.     From 
equations   (6)  and  (7)  we  can  compute 


will/)  ■     ^x     ay 
J(-s   ;       ax'  S  y' 


3x 


^y 


3y!     as' 


We  find  that,  at   the  origin  (i.e.   x'  =  y«  =  z»  =  0),   the  above  Jacobian 
is  equal  to  one.     Using  this  result    in  the  preceding  expression  for 

J(— )  we  obtain 
P 

J(^)    =  J(sj   . 

P 

With  the  x-,y-,z-axes  chosen  as  above  we  may  calculate,  by  using  equations 
(h)   in  equation  (5),    th-6  above  Jacobian: 

a?  viT  ;       aT7  vii  ; 


j(2l)  =  j(l)  -Jl  < 


TI 


+    Z 


9   <2^ 

ax1  **5v 


it 


ar 


^       ay'    ^     ax»%ij." 


-1 


■z     -  -     *z         *z 
where  the  partial  derivatives  are  to  be  evaluated  at  the  origin.   The  com- 
ponents of  the  vector  T  and  their  partial  derivatives  evaluated  at  the 
origin  are: 


T     - 
z 


1, 


x 


I 
.1-1 


1  - 


1/2 

•       =  cos  o( 


-  n 


V\ 


a1"1  VU\ 


3i 


::         n 


i-l 


(l-~)-2ai 


S»%_      i-l    Vi 


a  y 


-  -II 


"? 


C9) 


3% 

a  y 
91*, 

3  * 


i-l 


ri*  ®\=-  ni-i  !i!i 


(1  -  -ft-)  -  2b  A     , 


-  2a  An 


1-1  *1     +  n2(i-l) 


.:•  - 


-: 


i-i 


i  J 


i 


cos  « 


%, 


„     ,   i-l       yl      .      2(1-1)     ylY 
-2b  An  -g-    +  n 


? 


where  we  have  used: 


A 


.«  J  -(-l)1  {l  -  n^"1'  H   . 


2  -,\L/2 


>?-" 


(10) 


a      *  cos1^   +  cos  Of.  ; 


D2     =  s2  +  v2  +  z2 
1         T.       yl         1 


2^ 
cos    j 


and  the  definition  of  n.     With  the  aid  of  equations    (9)  the  JacoMan  in 


s 


(8)  'becomes 


P 


1  + 


COSCu 


i-l/      4      ^(i-D^^N 


1  +  —  + 


4        2 
D/cos    CC 
1  i 


-2(a+"b)A 


ax?  +  "by.2       2(i-l) 
-2(    ^  0       X)   r'       „ 


D- 


2 

cos  c< 


2 

cos  o(  1 


n2^/    «?     +  a^-^jUf+yM _  ^« 


d; 


-2(a+lD)(2~-A 

Dl  cos^O^        y  \^l 


-2(— 5 )     — 2~  K-Aj 

D-j_  cos  a.    >>    1      J 


+4at>£l  + 


— _ -p JA     +  2( ^ hi = 11  _- —  a 

D^  cos  «i  /  D^  v        D^cos^oc  /V>1       /J 


.-1 


This   equation  may  "be  expressed  in  terms  of  geometric  properties  of  the  sur- 
face and  distances  which  are  independent  of  the  choice  of  coordinates. 
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These  quantities  are 


(12) 


G =  (a  +  d)  ;        G     =  hah 

IU  g 

G^=  2(axf  *  oyr      ,       1 


hyf 


) 


sin2V 


V    2     v  2 

tx!  *  ay,  i 

2Gn  -  gt1  =  2(    ^     -    X  )  — 3y- 
m         U  D2  sin2r 

and  the  distances   involved  are  D.  .already  given,   and 


(13) 


D2  =  x2  +  y2  +  z2  = 


cos  c*. 


G  and  G  axe  respectively  the  mean  curvature  and  the  Gaussian  curvature 

of  the  surface  at  the  point  involved;  G^jis  the  curvature  of  the  surface 
in  a  plane  containing  the  incident  ray  and  the  normal  to  the  surface  at 
this  point;  D^  is  the  distance  along  the  incident  ray  from  the  source  to 
the  surface,  and  D  is  the  distance  along  a  reflected  or  transmitted  ray 
from  the  surface  to  any  point  (x,y,z)  on  this  ray.   In  terms  of  these 
invariants  equation  (11)  becomes 

f    r  i-1 

J(*r)  =U  +  SjV^1  *  ££lT-)  "   (2Gm  *  Gu^^K^cosrV  cos  ^)| 


J  n1"1,,        cos2Y\ 

'I1  *  ^  *  Z^ 


(HO 


2 
+  D 


a2(-  ^    /cos2y\       ,„  „      .      2r>    /cos^r\/i-l        y.  ^n1'1 

2 \ 2~7  "  (2&m  +  GlLtan  °    V 2— An       cosr-coscc)  — 

u     3^  cos  c^  cos  ^  ^i 


G     sec    o<.(n"""  cos  tf*  +  cos  o( .  ) 

g  i  i 


-1 


ITow  with  the  origin  at  any  point  in  space,   let  the  source  be  at 
(z1,y1,si),   let  a  point   on  the  surface  "be     (x1  ,y '  ,s' ),and  finally  let  the 
observation  point  be   (x,y,z)   in  the  new  coordinates.      Equation   (lU)   is 
unchanged  in  these  coordinates  provided  that 


(15) 


D^  =    (x1-  x»)2  +   {Vl~Vl)  *  (zx-  z')2 


D2  =  (x-x«)2  +  (y-y')2  +  (z-z«)2 


(z-z»)2 


2 
cos    ex. 


If  equation  (1U)  is  used  in  equation  (l)  the  reflected  or 
transmitted  field  amplitude  is  given  by: 


(16) 


E(x,y,z)  =  E^x'.y'.z')  J1/2(&mU'  ,y '  ),G  Jx'  ,y « ),(i]1  (x>  ,y» ) ,  ^,  Y^.d) 
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since  we  can  express  x  and  y  in  terms  of  x',y'  and  s  "by  means  of  equations 

(h).  In  equation  (l6)  J  is  a  functional  symbol  for  the  expression  on  the 
right  side  of  equation  (lh)  with  the  changes  noted  in  (15). 


h.   Discussion  of  the  Solution 

It  is  of  interest  to  investigate  the  "level  surfaces ,B  that  is 
surfaces  on  which  a  reflected  or  transmitted  field  component  has  a  con- 
stant value.  To  find  such  surfaces  we  require  that 

E(x,y,z)  =  3^  =  const, 
in  equation  (l6).  We  then  solve  the  resulting  equation  for  D  and  ohtain: 

f   /2?(x',y',2')\]  frrfl  +  r/H2    \  E?(x»fyiiSt)  11/2]  " 
(17)  D(x9y,z) 

') 

2 

where  H  and  K  are  the  coefficients  of  D  and  D  respectively  in  equation 

(lU).  Equations  (U)  and  (17)  are  equations  for  the  level  surfaces  in  terms 
of  the  parameters  x'  and  y'.  For  each  value  of  E„  there  are  two  such  sur- 
faces as  the  double  sign  in  equation  (17)  indicates. 

Those  level  surfaces  on  which  E  is  infinite  are  called  caustics 
and  are  of  special  importance.  Their  equations  are  easily  obtained  fcy 
letting  E#-*  00  in  the  equations  for  the  level  surfaces.  The  distance 
from  a  caustic  to  the  surface  along  a  reflected  or  transmitted  ray  is 
useful.  This  distance  is  given  hy  (17)  with  E  =  00  and  is: 


1  - 1— W. — ))\ W  K  -  ^      4      KJ 


(IS) 


D.|-|  iV(|)2-K   j 


thus  each  reflected  or  transmitted  ray  intersects  two  caustics. 

The  two  points  along  a  reflected  or  transmitted  ray  at  which 
E(x,y,z)  =  00  are  called  conjugate  points.   If  the  caustic  surfaces 
intersect  then  the  conjugate  points  on  the  rays  through  the  intersection 
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coincide,   and  a  point   image  is  formed;      the  image  is  called  real   if  D  >  0 
and  virtual  if  D  <  0.      This  occurs   if  the  second  term  in  (IS)   is   zero 


that   is 


IT  -  Iff 


_     n 


2(i-l) 


2  Y~ 

cos       » 

2    . 
cos    oC, 


-  1 


(19) 


+2 


i-1 


h 


2Gm  -  Gntan  c<i 


2^ 

cos     d 

2 
cos  c<. 


-  1 


i-1 


(n       cos  1+  cos  c*.  ) 


+   (2G  -K}-,,tan2<*  )  _1«S  sec2c*        fni"1Cos  *+  cos<*.  )2  =  0 
_     a     -.J.  x  g  x      x  i 

and  corresponds  to  the  following  possibilities: 

For  reflection,  o<  =  <f( these  results  are  also  given  in  ref.   3): 


l)  f=    tt/2  ;      then  from  eq.(lS)  D  =  -  D.,    and  the  virtual 
image  is   the   source. 
2)  a  =  "b  =  0;      all  possible  angles  for  V  and  D  =  -  D 

This   is   the   case  of  reflection  from  a  plana 
surface. 


3)     a  >  h  >  0 
a  <  h  <  0 


,   cos  V"  =   ■yfbfa  ,     9     =  tt/2 


U)     h  I  a  >  0 

h  <  a  <  0 

?or  transmission: 


,  cos  r  =  ya/D     ,0    =    0 


1)  n  =  1;      then  from  eq.    (IS)     D  -  -  D.    and  the  virtual  image 
is   the  source. 
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2)  Y*=  0,   c*2  -  tt;      then  Gm  =  Gg  and  the  point  of  trancnission 

Is  an  umbilical  point.      From  equation  (IS)  the  distance 
is  given  "by 

*  =  -     2_+  2a(n-l) 
D 

In  all  of  the  above  cases  the  distance  to  the  image  is  expressed  "by  the 
general  formula 

(20)  i    --    %    '[l-C-U1^1/2 

where  the  sign  to  he  chosen  is  that  of  G  .  If  the  source  is  at  an  infin- 
ite  distance  from  the  reflecting  or  transmitting  surface,  the  distance  of 
a- resulting  point  image  (vhich  is  called  a  focal  point)  from  the  surface 
is  called  the  focal  length,  f ,  and  is  given  "by  equation  (20)  with  Dt=  oa  . 


/2 


(21)  7  ■  *  t1  "  (-1)i  a]  GgX 

From  equations  (20)  and  (21)  we  obtain  the  combined  lens  and  mirror  law: 

(22)  -  +  —  =  I 

As  another  application,  we  specialize  the  reflecting  surface 

to  a  sphere.  Then  equation  (lh)  yields  the  geometrical  factor  obtained 

10 
by  van  der  Pol  and  Bremmer  for  the  reflection  of  an  electromagnetic  wave 

from  the  earth.  V/hen  the  Presnel  formulae  and  the  phase  are  taken  into 

account,  their  complete  result  (for  wavelengths  small  compared  to  the 

earth's  radius)  is  obtained. 
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APPETDIX 


Derivation  of  Eouation  for  T 


To  derive  equation  (2)  of  the  test  proper  we  use  the  unit 

vectors  and  angles  shown  in  Fig.  1.   From  the  definitions  of  these 

vectors  and  angles  ve  require: 

->2   -» 2    i-?2 

I  =  r  =  TT  =  1 

/r\  I*  •  H  =  cos  >^ 

5  •  N  =  cos  <* 
By  the  law  of  reflection  V"=o',   or: 

(7)  3'T   •  1?  »    "i   •   H     ; 

and  from  Snell's  law  n  =     sin  ^2/   sin  V*    ,   or: 

T?  •  H  =     -  |l  -  n^  [l  -  (I  •  N)2J  J 

where  we  always  take  the  positive  square  root  in  Eq.  (3).  Since  the  re- 
flected or  transmitted  ray  must  lie  in  the  plane  of  the  normal  to  the 
surface  and  the  incident  ray  (by  the  law  of  reflection  and  Snell's  law) 

we  have 

—  i->     — >     — * 

00  T?  =  Ajl  +  3^^  H  . 

Using  equation  (4)  in  equations  (1)  and  (2)  we  find  the  two  solutions 


The  second  set  of  coefficients  reproduces  the  incident  ray  and  so  nust  be 

1-** 

discarded.  Then  the  most  general  T  vector  is  given  "by: 

1-5"     — *■     .-*■   ->.  -* 
(5)  T  =  -  I  +  2(1  •  ST)  N. 


ih 


If  equation  00   is  used,  in  equations   (l)  and(3)  we  get   the  soluti 


ons 


A2  =  -  n  and  A    =  n 

V2 
J32  =  n(f  •  f)  -  {l  -  n2[l  -(?•  tf)2]}  B2  =  -ntf  •   N)-^l-n2[l-(?.S)]] 

The  second,  set  of  coefficients  yield  a  transmitted  ray  which  lies  on  the 
same  side  of  the  normal  as  the  incident  ray.     However,   this  is  in  contra- 
diction to  the  law  of  refraction  and  so   this  solution  must  he  abandoned. 
Using  the  first  set  of  coefficients  gives  the  general     x  vector  as: 

(6)         ¥--«?*(ir-t)-{i-a2^-(?.i)g}    )i 

Equations  (5)  and  (6)  can  now  he  written  in  the  comhined  form: 

(7)    S  «  -(n1"1)  T  ♦  (»"(?  •  SM-ltfl-n2'1-1^.?)2]}172)  t 


which  is  the  desired  equation.  This  can  also  he  expressed  by  the  simple 
expression: 

(8)  ^T     =   (-ni_1)  t  +  (n1"1  cos^*  cos  oC  )  f 


1/2 
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